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INTRODUCTION 


The interaction of atomic particles with surfaces is of both scientific and 
technological interest. Past work emphasizes the measurement of high-energy sputtering 
yields. See Behrisch, R., editor, Sputtering by Particle Bombardment I: Physical 
Sputtering of Single-Element Solids, (Springer-Verlag, Berlin, 1981) for a survey of 
published sputtering yield data. Very little work utilized low-energy beams for which 
chemical and electronic effects can be important. Even less work has been carried out 
using well-defined low-energy projectiles. The use of low-energy, reactive projectiles 
permits one to investigate surface processes that have not been well characterized. As the 
energy of the projectile decreases, the collisional cascades and spikes, that are common 
for high-energy projectiles, become less important, and chemical and electronic effects 
can play a significant role (Sigmund, 1981). 

Aspects of particle-surface interactions are of concern in several areas of 
technology. For example, the erosion, desorption, and glow of surfaces of spacecraft in 
orbit are important in the arena of space technology. The materials studied under this 
contract are of possible use on the exterior portions of the power generation system of 
Space Station Freedom. Under the original designs. Space Station Freedom's power 
generation system would generate potential differences on the surface as high as 200 
volts. Ions in the plasma that often surround orbiting vehicles would be accelerated by 
these potentials leading to bombardment and erosion of the exposed surfaces. The major 
constituent of the atmosphere, approximately 90%, in the low earth orbit region is atomic 
oxygen. Since atomic oxygen is extremely reactive with most materials, chemical effects 
can arise in addition to the physical sputtering caused by the acceleration of the oxygen 
ions. Furthermore, the incident oxygen ions can remain embedded in the exposed 
surfaces, altering the chemical composition of the surfaces. Since the effective binding 
energy of a chemically altered surface can be quite different from that of the pure 
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substrate, the sputtering yield of a chemically altered surface is usually different also. 
The low-energy 0 + sputtering yield measurements, reported here, will help quantify the 
erosion rates for materials exposed to the low-earth orbit environment. 

These measurements are of technological importance in another respect. In most 
surface analysis techniques, a surface is bombarded with ions, electrons or photons. 
Information concerning the structure of the surface and near-surface bulk, abundance of 
impurities and defects, as well as other surface properties are obtained either from the 
desorbed species or from the scattered projectiles. Because of their low penetration 
depth, low-energy ions provide an advantage over other techniques because they provide 
information that is more indicative of conditions on the surface rather than integrated 
effects arising from deeper in the bulk. A better understanding of the microscopic 
processes involved in these interactions is not only of basic scientific interest, but will 
also aid the scientific community by increasing the accuracy and usefulness of these 
surface analysis techniques. 



EXPERIMENTAL METHODOLOGY 


General 

The erosion rate experiments described in this report were conducted in an ultra- 
high vacuum chamber with a base pressure of 1 x 10"9 Torr (1 x 10'8 Torr when running 
an ion beam). In these experiments, a low-energy ion beam was incident on a thin-film 
sample that covered a quartz-crystal microbalance. Bombardment of a surface with a 
low-energy reactive ion beam can simultaneously remove and add material to the surface 
of the substrate - sputtering of the substrate and incorporation of projectile particles into 
the substrate. Since the samples were thin films deposited on the surface of a 6 MHz AT- 
cut quartz crystal, the net mass change induced per incident ion was obtained by 
measuring the induced frequency shift of the quartz crystal. The current on a nickel grid 
was simultaneously measured to determine the number of ions incident on the sample. 
These measurements provided the information necessary to calculate the sputtering yield. 
Details concerning the apparatus and experimental procedure are discussed below. 

The Ion Source and Accelerator 

The ion beam is produced by a modified commercial Colutron ion source. The 
oven, where the gas is ionized and a plasma is maintained, is essentially a hot cathode 
type ion source (Figure 1). A quartz tube serves as the containment walls of the oven. At 
the back end of the oven, the gas inlet tube introduces the beam gas at pressures that are 
generally in the hundreds of mTorr region. The front of the oven is closed except for the 
0.020 inch anode hole. This arrangement allows differential pumping between the inside 
of the oven and the rest of the ion source, where operating pressures are usually in the 
flTorr region. 
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To create an ion beam, the gas atoms or molecules must first be ionized and even 
disassociated in cases where the introduced gas is a molecular species and atomic ions are 
desired, e.g. 0 + from O 2 . A beam of electrons ionizes, and in some cases dissociates, the 
gas. The filament is heated by a current to a sufficient temperature to "boil" electrons 
from its surface. These electrons are then accelerated toward the anode by a 40- to 100- 
volt potential difference maintained between the anode and the filament. The positive 
ions that result from electron/atom collisions and ion/atom collisions form a plasma near 
the anode. Positive ions near the anode hole are then extracted by a negative 2-keV 
potential difference maintained between the extraction plate (see figure 2) and the anode. 

The very corrosive, high temperature environment created in the oven when used 
to create an oxygen plasma necessitated the replacement of all the metal parts in the oven 
with more robust materials under these conditions. For the ion source, Colutron supplies 
tantalum anodes, alloy filament rods, and tungsten filaments. In a high temperature, 
oxygen- rich environment, the anodes oxidize and deform, the filament posts oxidize and 
coat the inside of the oven with metal oxides, and the tungsten filaments bum. The 
oxidation of the tungsten quenches the plasma and stops the discharge. To address these 
problems we replaced the anode and filament posts with similarly shaped platinum 
pieces. The platinum pieces are much less reactive with the oxygen and experience very 
little degradation even under these adverse conditions. Since the platinum filament posts 
become soft and move when heated, we designed a quartz ring with the appropriate holes 
to support the posts. The tungsten filament was replaced by an iridium ribbon coated 
with thorium oxide to reduce its work function. The iridium filaments last more than an 
order of magnitude longer than the tungsten in the oxygen environment. 

The positive voltage applied to the anode determines the final beam energy (beam 
energy = anode voltage). The extraction plate is held 2 kV beneath the anode voltage to 
maintain a high extraction efficiency and to minimize the effects of space charge. The 
power supplies that control the beamline are then "floated" at the same potential as the 
extraction plate (Figure 2). All of the beam optics are referenced from this level and are 
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surrounded by a wire grid that is floated at the same potential in order to create a field- 
free region for the beam. Even though the ions are traveling with velocities defined by 
the 2-kV accelerating potential, they only have a potential energy, with respect to ground, 
that is equal to the voltage on the anode. We are therefore able to decrease the kinetic 
energy of the beam to the desired amount by placing two nickel grids (90 percent 
transmission, 70 lines per inch) in the path of the beam (Figure 4). The grids are 
perpendicular to the beam axis and parallel to each other. The first grid is held at the 
accelerating potential, while the second grid and the sample are grounded. The beam 
essentially doesn't "see" the first grid nor can it "see" the second grid until after it passes 
through the first grid. After passing through the first grid, the ions are decelerated. The 
kinetic energy of the ions is thus reduced to an energy determined by the voltage on the 
anode. 

The extraction plate is mounted on the first section of an Einzel lens located 
immediately after the anode. The Einzel lens focuses the beam which is diverging as it 
emerges from the anode. As depicted in figures 2 and 3, the beam then goes through a set 
of vertical deflection plates for small steering corrections and a Wien filter for 
mass/velocity selection of a given charge state. The Wien filter allows selection of either 
an 0 + or an C> 2 + beam, for example, from the mixture that is extracted from the anode. 

The primary purpose of the four sets of deflection plates that are encountered after 
the Wien filter is to remove high-energy neutrally charged particles from the beam. The 
neutral atoms or molecules arise from collisions during or immediately after extraction 
from the anode. The first set of plates is oriented in the vertical direction and serves to 
bend the beam to the right. The second and third sets of plates are positioned horizontally 
and are used to make slight corrections to the beam's vertical trajectory. The fourth set of 
plates, oriented in the vertical direction, bends the beam back to the left along a direction 
that is parallel to the original beam axis, but is offset 0.5 inch to the right. A copper 
block (l"xl"xl") with a 3/16" hole, angled at 1.3 degrees with respect to the original 
beam axis, is located between the second and third sets of plates. The copper block is 
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mounted on a plate that separates the high vacuum ion source from the ultra-high vacuum 
sample chamber. The small hole in the copper block allows for differential pumping by 
the two 200 liter/sec ion pumps located on the ultra-high vacuum side. 



Figure 2. Schematic diagram of high vacumm portion of beamline. 
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A second Einzel lens, located after the fourth set of deflection plates, focuses the 
beam through the beam tube and onto the sample. The beam tube is floated at the 
accelerating potential and serves to shield the beam from the walls of the bellows, which 


Figure 4. Top view of final portion of beamline and sample chamber. 
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are grounded. Figure 4 shows that we have placed a permanent magnet between the last 
deceleration grid and the sample. When the ion beam strikes the grids, secondary 
electrons can be ejected from the grid. The electrons that are created at the first grid are 
then accelerated towards the target by the same potential difference that decelerates the 
ion beam. The permanent magnet diverts the electrons away from the target while 
bending the positive ion beam only slightly in the opposite direction. 

Using the Ouartz-crvstal microbalance for Mass Loss Measurements 
To measure erosion rates or sputtering yields, one needs to count the average 
number of particles emitted from the sample surface when a projectile ion strikes the 
sample at a given energy. This task is usually accomplished by measuring the total mass 
change of the sample when it is struck by a large known number of ions. From these data 
the total number of emitted atoms (or molecules) per incident particle is calculated. 
Because of its high sensitivity, we used a quartz-crystal microbalance (QCM) to measure 
the mass change. The samples were thin films sputtered or evaporatively deposited on 
the surface of the quartz crystals. 


Gold electrodes 



6 MHz AT-cut 
plano-convex 
quartz crystal 


Figure 5. Side view of the quartz crystal. 



OCM Theory and Operation 

A quartz-crystal microbalance is comprised of a thin quartz wafer sandwiched 
between two metal electrodes, an oscillator circuit, and a circuit for measuring the 
resonant acoustic frequency of the crystal/thin-film system (see Figure 5). The oscillator 
circuit is utilized to apply a radio-frequency voltage across the electrodes which drives an 
acoustic mode in the piezoelectric quartz. This resonant acoustic mode is extremely 
sensitive to any mass loss or gain from the surface of the crystal. Sauerbrey was the first 
to methodically investigate the possibility of using a quartz crystal in the configuration 
described above to measure small mass changes (Ullevig and Evans, 1980). He noted 
that the mass change was proportional to the observed frequency shift; specifically, the 
magnitude of the mass change occurring on the surface of a quartz-crystal microbalance 
can be determined with the aid of Sauerbrey's relationship 


A/ = - 


2/qA/m 


(i) 


where Af is the observed frequency shift, /„ is the parent frequency of the quartz crystal, 

Am is the mass change, A e is the active area of the crystal (i.e.. the area of the smallest 
electrode), p q is the density of quartz (2.648 g cm 3 ), p q is the shear modulus of the 

crystal (2.974 * 10 1 1 dynes cm 2 for an AT-cut quartz crystal (Ward and Delawski, 
1991)). 

The constants in equation (1) can be grouped to obtain the relationship 

A/ = ~ C 7 ,Am (2) 

where C lola , is often called the sensitivity constant for the crystal. It should be noted that 
equation (2) is valid only if the induced frequency shift is small in comparison to the 
initial frequency of the crystal and the mass is added or removed uniformly across the 
active area of the crystal (Ullevig et ai, 1982). 
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If the mass change occurs over an area that is smaller than the active area of the 
crystal, the magnitude of C will depend on the area distribution over which the mass 
change occurs and the sensitivity function, S(r), which describes the incremental 
frequency shift, df, induced by an incremental mass change, dm, at a given distance from 
the center of the active area. Stated differently, the sensitivity constant is merely a 
weighted summation of the sensitivity function across the area where the mass change 
occurs. Therefore when the mass change occurs uniformly across the total active area of 
the crystal, the sensitivity constant is 

C ,o,at =27r V r S{r)dr where S(r) = ^j^-. (3), (4) 

J0 dm 

The subscript "e" indicates that the entire electrode area is utilized. If the same amount of 
mass, Am, were removed uniformly from a spot that is smaller than the active area, then a 
different frequency Af' would arise and the integration should be carried out over the spot 
area. Similar to equation (2), the new frequency would have the form: 

— C .Am 

A r=— f 2 where C spot = 2 k)~ rS(r)dr. (5), (6) 

™spot 

From an intuitive standpoint, it is not surprising that S(r) is largest at the center of 
the crystal and decreases monotonically to the edge of the active area since the crystal is 
vibrating in a thickness shear mode. Several groups have measured the local vibration 
amplitudes of the quartz crystal - Sauerbrey (1964) used the modulation of a light beam, 
Wimmer et al (1984) made use of the speckle effect, and Mecea (1988) used a small 
rubber tipped wire to measure the vibration energy dissipation which he then converted to 
vibration amplitude. These measurements revealed that for an AT-cut quartz crystal, the 
amplitude of vibration is largest at the geometric center of the crystal's active area and a 
direct correlation can be made between the amplitude of the acoustic wave at a given 
concentric ring and the magnitude of the sensitivity function for that ring. 

An Inficon XTC quartz crystal controller was used to drive the acoustic mode of 
the crystal and to measure the frequency shift induced when the ion beam removed mass 
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from the surface of the crystal. The controller automatically reports a mass change 
which corresponds to the frequency shift that would be induced if the mass were removed 
uniformly across the entire active area. However for this experiment, we chose to 
uniformly remove mass from a spot smaller than the electrode area in order to concentrate 
the small currents obtainable with our low-energy ion beam on the most sensitive area of 
the crystal. It was therefore necessary to relate the observed frequency shift, Af\ to the 
frequency shift, Af, that would have been induced by the same change in mass. Am, if the 
mass were removed from the entire active area. To obtain this relationship, we divide 
equation (1) by equation (5) after replacing the sensitivity constants with the appropriate 
integrals from equations (3) and (6) to obtain 



fact, Ullevig et al. (1982) showed that the integrated sensitivity constant is constant for a 
given crystal cut, but the sensitivity function will vary depending on the physical nature 
of the thin film on the crystal's surface. Furthermore, we only have the ability to measure 
a quantity, which we call F(r), which is a convolution of the true sensitivity function with 
the finite area distribution from which the mass is removed. Measurement of the true 


sensitivity function would require removing mass from an infinitely small area, which is 
of course impossible with an ion beam of finite dimensions, while simultaneously 


measuring the induced frequency shift of the crystal. In fact, we determine F(r) by 
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measuring a quantity that is proportional to df/dm. Fortunately the proportionality 
constant is not needed since F(r ) is used in both the denominator and the numerator of the 
modified version of equation (7). The procedure for measuring F(r) will be discussed in 
a separate section. In the following section, we will show that the ratio of Af to Af' given 
in equation (7) can also be computed using F( r). 



Figure 6. Active area of crystal - radius r e ; spot radius r S pot- 

First, consider Figure 6 where the larger circle with radius r e represents the active 
area of the crystal while the smaller area with radius r S pot represents the area of a hole in 
the mask that covers the crystal limiting the area of mass removal. As noted earlier, F(r) 
is actually a convolution of the true sensitivity function with the area from which the 
mass is removed, or 

F(r„) = \“\frSVX)drde (9) 

where SR is the distance from the center of the active area to an infinitesimal area element 
of mass removal within the spot. As shown in Figure 6, it can be expressed as 

SR = o + 2rr 0 cos 6 + r 2 . (10) 

We noted earlier that F(r) is a convolution of the true sensitivity function with the 
finite area distribution from which the mass is removed. Closer inspection reveals that 
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the sensitivity constant can be computed using F(r) if the integral is normalized by the 
area of the spot, i.e. 


2ftj o ' rS(r)dr 


2/rj*^ rF(r 0 )dr 


nr. 


spot 


( 11 ) 


Furthermore, Equation (10) shows that a special case exists when the spot is centered on 
the active area, i.e. when ro is zero. 


r 2 k r r inn . e rspot 

F(0) = j \ Q r o S(3i)dr o d0 = 2n) o rS(r)dr (12) 

Substituting equations (1 1) and (12) into equation (7) brings us to our objective - relating 
Af to Af using F( r) instead of S( r). 


2 nj ' rF{r)dr 

Af (^Xgj ifcrFi'Vr 

A f m r 2 F( 0) K 

^ r spot 

It is now possible to find the true mass change from the reported mass change 
when the mass is uniformly removed from a spot that is smaller than the active area of the 
crystal using the relationship 


Am = Am 


reported 


4f 

A/'' 


(14) 


Measurement of the Convoluted Form of the Sensitivity Function 
To measure F( r), the QCM is first positioned so that the beam strikes the center of 
the active area while the area of mass removal is limited to the size of the hole in the 
mask. To insure uniform mass removal across the small spot area, the ion beam is 
rastered over an area slightly larger than the hole in the mask which defines the spot area. 
As described earlier, the number of ions striking the sample is recorded along with the 
reported mass change, Am ' from the QCM controller While the beam and the mask 

remain stationary, the QCM is then moved so that the center of the small spot is 
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positioned a small radial distance from the center of the active area. Again the number of 
incident ions is recorded along with the reported mass change which should be smaller 
than the first since more mass is removed in this case from areas having smaller S(r). 
This procedure is repeated - effectively moving the area of mass removal across the 
diameter of the active area. The magnitude of Am 7(# incident ions) at a distance r from 
the center of the active area is equal to F(r) which is proportional to df/dm since the 
reported mass change is actually an interpretation, by the QCM controller, of the induced 
frequency shift. Furthermore, the normalization to the number of incident ions is 
equivalent to a mass normalization since an incident ion on the average always removes 
the same amount of mass regardless of the incident location. A plot of the results from 
this procedure produces a curve like the one in Figure 7 which is proportional to the true 
sensitivity function of the crystal S(r). 


T3 

o 

u 


3 

Cfl 

a 

<u 


<u 




a 

E 


u 

© 


Z, 



Figure 7. The convoluted sensitivity function vs. the crystal radius. 
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Sample Preparation 

The gold electrodes on each side of the quartz are approximately 2500 A thick 
and were evaporatively deposited by Inficon on a thin layer of chromium which bonds to 
the quartz better than gold (Roekel, LV., personal communication, February 19, 1992). 
The NASA Lewis Research Center prepared the aluminum and aluminum oxide samples 
by sputter deposition on the larger gold electrode as shown in Figure 5. The thickness of 
the aluminum oxide film was 1500 A and that of the aluminum film was 2200 A. 
(Vaughn, J., personal communication, February 19, 1992). 

Prior to the erosion rate measurements, each sample was sputter-cleaned by a 2- 
keV argon-ion beam. A minimum of 40 A of material was removed; and the argon 
bombardment was continued until the measured erosion rate of the sample became 
constant. 


Data Acquisition 

We determined the number of ions striking the sample by measuring the current in 
a grid located behind the 2-mm hole (see Figure 4). Prior to the erosion-rate experiments, 
we determined an energy-dependent conversion factor that relates the current incident on 
the sample surface to the current measured on the grid. (Corrections were made for 
secondary-electron emission.) An Inficon thin film rate monitor was used to drive the 
quartz crystal and extract the frequency shift needed for the determination of the mass 
change. A Macintosh Ilci computer equipped with a National Instruments GPIB interface 
card was used to control the experiment, record the data, and perform the necessary 
calculations needed to extract the sputtering yield from the data. 

To obtain the sputtering yield from the raw data, we first performed linear 
regression on the frequency data which are a function of the total ion dose. We then 
compared each measured frequency to the value predicted by the equation from the linear 
regression procedure. Measured frequencies that were more than 3 standard deviations 
from the predicted value were discarded and linear regression was performed again on the 
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remaining points. The slope obtained from the second linear regression procedure was 
then multiplied by the appropriate conversion factors to transform the frequency shift 
(corresponding to a specific mass change) to the corresponding number of sample atoms. 
That result was multiplied by another conversion factor to change the denominator to the 
number of incident ions, giving an uncorrected sputtering yield. This value was corrected 
by the factor of equation (14) to give our measured value of the sputtering yield. 



EXPERIMENTAL RESULTS 


The thicknesses of the films supplied by the NASA Lewis Research Center are shown in 
Table 1. 


NASA Samples 

Table 1 . Approximate film thicknesses of the NASA samples. 


NASA 

Samples 

Approximate 
Film Thicknesses 

0 

(A) 

Cu 

800 

SiO x 

1800 

AI2O3 

1500 

A1 

2200 


Copper Sputtering Yield Results 

For the copper sample, we measured the sputtering yields for low-energy Oj + , 
0 + , N2 + , and N + bombardment. The results are presented in Figure 8 and Table 2. The 
solid curves shown are second-order polynomial fits. Several low-energy sputtering 
measurements were made for N2 + on copper to determine the sputtering threshold. As 
can be seen in Table 2, the lowest energy at which a mass loss of the copper was 
measured was 60 eV. No measurable mass loss was observed for a 55 eV N2 + beam. 

For oxygen bombardment of copper, a definite dependence of the measured 
sputtering yield on bombardment time was observed for both oxygen ion species. This 
time dependence is probably due to both an initial mass gain caused by the reactive 
incident particle and the true difference between the sputtering yields of the oxides 
formed on the surface layer and the sputtering yields of the pure metal. To quantitatively 
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address this issue, we devised a cleaning procedure that is discussed below. Preliminary 
experiments with the copper indicated that the time dependence of the C> 2 + and 0+ 
sputtering yields not only depends on bombardment time, but also on whether the copper 
surface was cleaned with an argon beam before oxygen bombardment. The copper 
samples were sputtered through to the gold substrate before more complete measurements 
could be made. 



Figure 8. Sputtering yield results as a function of energy for copper. 
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Table 2. Sputtering yield ( average number of ejected atoms per incident ion) 
measurements forN2 + , N*, Of*, and O* on copper. 


Projectile 

Energy 

(eV) 

n 2 + 

N+ 

Of* 

0+ 

50.0 


0.037 



60.0 

0.022 




65.0 

0.027 




75.0 

0.045 




100.0 

0.14 

0.15 

0.015 

0.052 

150.0 

0.24 

0.25 

0.092 

0.11 

200.0 

0.32 

0.33 

0.17 

0.16 

400.0 


Billillii 

0.57 



Silicon Oxide Sputtering Yield Results 

For the silicon-oxide sample, we measured the sputtering yields for low-energy 
C> 2 + , and 0 + . The results are plotted in Figure 9 and listed in Table 3. 

Because the exact compositions of the SiO x samples are unknown, mass losses 
are reported instead of sputtering yields. Beam focusing and distortion was a problem for 
all measurements due to charging effects of the non-conductive sample. The mass losses 
are small for all bombardment energies studied. For example, a mass loss of 1.99 x 10 22 
g for one 2 keV C> 2 + ion corresponds to a sputtering yield of only 0.42 if the sample was 
pure SiC >2 or 0.57 for pure SiO. 


Table 3. Mass loss measurements for Of* and O* on SiOx. 


Projectile Energy 
(eV) 

0 2 + (x lO- 22 g) 

0+ (X 10-22 g) 

250.0 

0.095 


500.0 

0.38 

0.35±0.05 

1000.0 

1.01 

1.02+0.07 

1500.0 

1.41 

1.35+0.26 

2000.0 

1.99 































Mass Loss (x 10 22 g) 
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Aluminum Oxide Sputtering Yield Results 
For the aluminum-oxide sample, we measured the sputtering yields for low- 
energy C> 2 + , 0 + , N 2 + , and N + bombardment. The results are plotted in Figure 10 and 
listed in Table 4. 



Sputtering Yield 
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Energy of Incident Ion (eV) 

Figure 10. Sputtering yields as a function of energy for aluminum oxide. 
Table 4. Sputtering yield measurements on aluminum oxide sample. 


Projectile Energy 
(eV) 

Projectile 

Y measured 

(avg. # ejected atoms per 
incident ion) 

2000 

Ar + 

1.1 ±0.3 

400 

n 2 + 

0.41 ±0.10 

300 

n 2 + 

0.22 ± 0.06 

200 

n 2 + 

0.06 ± 0.02 

150 

n 2 + 

0.03 ±0.01 

400 

N+ 

0.33 ± 0.08 

300 

N+ 

0.19 ±0.05 

200 

N+ 

0.04 ±0.01 

400 

o 2 + 

0.29 ± 0.07 

300 

o 2 + 

0.17 ±0.04 

200 

o 2 + 

0.07 ± 0.02 

150 

o 2 + 

0.03 ±0.01 

100 

o 2 + 

0.01 ±0.004 

400 

o+ 

0.28 ± 0.07 

300 

o+ 

0.17 ±0.04 

200 

o+ 

0.07 ± 0.02 
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Based on what was learned from the Cu and SiO x samples, a two step cleaning 
procedure was developed and used for all AI 2 O 3 and A1 samples. First, the sample 
surface was bombarded with a 2-keV Ar + beam until more than forty monolayers of 
material were removed and the erosion rate reached a constant value. After the argon 
beam was used to clean the sample and to measure the sensitivity function of the crystal, 
the sample was unavoidably exposed to the residual gases in the sample chamber 
(pressure = 2xl0 * 9 Torr) while the desired ion beam was being tuned. During this time, 
atoms and molecules from the vacuum adsorb on the surface of the sample. These atoms 
and molecules, H 2 O, CO 2 etc., come from background gases in the sample chamber that 
stick to the surface during the time between the Ar + cleaning and the bombardment with 
the beam of interest. In the second stage of the cleaning process, the ion beam to be used 
for the actual sputtering yield measurement was first used to remove several monolayers 
of material from the sample while the sputtering rate was simultaneously measured. The 
second stage of cleaning was terminated when the measured sputtering yield became 
constant, indicating that the adsorbed species on the surface were removed. At the 
beginning of both cleaning processes, the mass change per incident ion is slightly larger 
than the final value. This difference occurs because the weakly bound surface adsorbates 
are much easier to sputter from the surface than are the aluminum-oxide molecules. The 
sputtering yields reported represent the constant values obtained after the surface 
impurities were removed. 


Aluminum Sputtering Yield Results 

The aluminum samples were cleaned by the procedure described above. During 
this procedure, it soon became apparent that the cleaning process was even more critical 
than it had been for the aluminum oxide. Again at the beginning of the cleaning process, 
the sputtering yield was larger than the final value. However, an intermediate stage was 
also observed in which the sputtering yield was smaller than the final constant value 
reached. Since the samples were exposed to air for several months before the experiment. 
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the top layers of the aluminum were at least partially composed of aluminum oxide. The 
measured results of the final constant sputtering yield reached are plotted in Figure 1 1 
and listed in Table 5. 

O — Ar + on clean A1 

• Ar + on oxygenated A1 



Energy of Incident Ion (eV) 

Figure 1 1 . Semilog plot of sputtering yield as function of energy for aluminum. 

The term "oxygenated aluminum" is used to describe the aluminum surface after 
bombardment with an oxygen beam. As will be discussed in greater detail below, the 
mass of the aluminum sample actually increases when the clean aluminum surface is 
bombarded with a low-energy ionic oxygen beam; and then the sample stops gaining 
mass and starts losing mass. The oxygen is obviously embedded into the surface and 
near-surface bulk of the initially pure aluminum. For oxygen bombardment of aluminum, 
we report the sputtering yield of the chemically altered surface obtained when the amount 
of mass loss per incident ion reaches a constant value. We are using the term 
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"oxygenated aluminum" to describe the surface of the aluminum sample after it has been 
exposed to a sufficiently large dose of oxygen that an equilibrium condition is reached 
and a net mass loss per incident oxygen ion is observed or no mass gain is observed. To 
investigate the effects of the chemically altered surface on the sputtering yield of a non- 
reactive bombarding species, a 2-keV Ar+ run was conducted on the oxygenated surface. 
The measured sputtering yield was much smaller than we had previously measured for 2- 
keV Ar + on the pure aluminum surface. 


Table 5. Sputtering yield measurements on aluminum sample. 


Projectile 

Energy 

(eV) 

Projectile 

Sample 

Y measured 

(avg. # ejected atoms 
per incident ion) 

Y literature 

(avg. # ejected atoms 
per incident ion) 

2000 

Ar+ 

A1 

1.6 ±0.4 


2000 

Ar+ 

oxygenated A1 

0.48 ±0.12 


1000 

Ar+ 

A1 

1.2 ±0.3 

2.00) 

400 

Ar+ 

A1 

0.69 ±0.17 

0.8(2) 

200 

Ar+ 

A1 

0.14 ±0.04 

0.35(2) 


WmsBM 

A1 

0.19 ±0.05 



■H 

A1 

0.13 ±0.03 



HSU 

A1 

0.07 ±0.02 


150 

n 2 + 

A1 

6.03 ±0.01 


400 

N+ 

A1 

0.18 ±0.05 


200 

o 2 + 

oxygenated A1 

0.03 ±0.01 


2000 

o+ 

oxygenated A1 

0.26 ± 0.07 


500 

0+ 

oxygenated A1 

0.06 ± 0.02 



1 (Oechsner 1973) 

2 (Laegreid & Wehner 1960) 


Dose - Dependent Oxygenation of Aluminum Results 
When a beam of 200-eV 02 + or 0 + ions was incident on a clean aluminum 
surface, a net mass gain per incident ion was initially observed. The oxygen was 
obviously being embedded into the surface and near-surface bulk of the aluminum. As 
the total dose increased, the magnitude of the mass gain per incident ion decreased. The 
mass change per incident ion finally reached a constant value after the samples were 
exposed to a very large total dose. For the runs having large enough beam currents 
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(making the mass decrease measurable), the final value was negative indicating a net 
erosion of material per incident ion. The open squares in Figure 12 correspond to a 200- 
eV 0 + run that exhibited this behavior. 

The aluminum sample was first cleaned with a 2-keV Ar + beam, as described 
previously. At the beginning of the run, the mass change per incident ion increased 
slightly from the first data point before it began the decrease toward the equilibrium 
condition. Implications concerning the equilibrium condition will be addressed below, 
but experimentally the equilibrium condition is defined as the point at which the mass 
change per incident ion becomes smaller than we can measure for a given run. We 
attribute the slightly smaller mass change per incident ion that occurred at the beginning 
of the run to the sputtering, i.e. negative mass change, of atoms or molecules that were 
weakly bound to the surface. As described above, these adsorbed atoms and molecules 
come from background gases in the sample chamber that stick to the surface during the 
time between the Ar + cleaning and the bombardment with the oxygen beam. This effect 
was observed for all runs of this type, and the magnitude of the initial increase in the 
mass change per incident ion observed corresponds directly to the amount of time that the 
adsorption of background gases occurred. 
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Cumulative Ion Dose (# of ions) 

Figure 12. Mass change per incident 200 eV 0 + ion on an aluminum surface. 


The run represented by the solid black squares in Figure 12 was performed on the 
same spot of aluminum that the first 200-eV 0 + run had oxygenated three days 
previously. The aluminum surface was not cleaned with argon prior to the second run. 
At the beginning of the second run, there was a positive mass change per incident ion that 
remained positive until considerably more oxygen was added to the sample. 
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Two additional runs were conducted in a similar fashion using 200-eV C> 2 + beams 
on two different aluminum samples. Unlike the procedure used for the 0 + runs, the 
sample surfaces were cleaned with argon before both runs. The major difference between 
the two 02 + runs was the beam current. The beam current for the run represented by the 
dark triangles in Figure 13 was roughly half that used for the run represented by the open 
triangles. As the figure shows, a positive mass increase per incident ion was observed 
for both runs. Furthermore, the mass changes per incident ion for the two runs are in 
good agreement for doses common to the two runs. 



200 eV Oj incident on clean A1 surface 
(14.3 ngrams of total mass gain before 
observable mass loss) 

Current = 40 nA 

200 eV Oj incident on clean A1 surface 
(1 1.2 ngrams of total mass gain before 
observable mass loss) 

Current = 75 nA 


. 0.0 


0.00e+0 


2.00e+l 5 


1 

4.00e+1 5 


1 

6.00e+1 5 


1 

8.00e+1 5 


Cumulative Ion Dose (# ions) 


Figure 13. Mass change per incident 200 eV (?2 + ion on an aluminum surface. 



DISCUSSION 


Dose Dependent Oxygenation of Aluminum 

When a 200-eV 0 + beam was incident on the clean aluminum surface, the mass of 
the aluminum sample was observed to increase for a surprisingly long period of time (see 
the data represented by the open squares in Figure 12). The positive mass change per 
incident ion, which arises because the incident 0 + ions are embedded into the surface or 
near surface bulk, is largest for small dose. After the aluminum was exposed to the 0 + 
beam for a long time, i.e. large dose, the mass change per incident ion finally decreased 
to a value that was smaller than we could measure. 

After the sample had been exposed to a large dose, an equilibrium condition was 
reached at which point the mass change per incident ion was constant. After the 
equilibrium condition was reached, the measured mass change per incident ion turned 
negative, it the beam current was sufficiently large. We believe that the equilibrium 
condition is an equilibrium between five competing processes: sputtering of the 
aluminum, incorporation of oxygen from the ion beam into the aluminum, sputtering of 
the incorporated oxygen, incorporation of oxygen from the background gas in the sample 
chamber, and sputtering of the adsorbed background gas atoms. The relative 
contributions of the five different processes are largely determined by the magnitude of 
the beam current. 

Mass gain due to the incorporation of oxygen ceases when the surface and near- 
surface bulk become saturated with oxygen. For the 200-eV 0 + bombardment the total 
mass increase of the aluminum sample before the mass change per incident ion reached 
zero was 13.8 nanograms. The total mass increase of a clean aluminum sample 
bombarded by a 200-eV C> 2 + beam was 14.3 nanograms, which is close to the value of 
13.8 nanograms measured for the 0 + run. (See Figure 13 and the accompanying 
description in the data section). This close agreement is also supported by other similar 
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runs. If we assume that no sputtering of the aluminum occurs during the period of mass 
increase, 14 nanograms corresponds to 11 oxygen atoms per aluminum atom in the 
surface layer. This is certainly not a valid assumption, but it sets a lower bound on the 
amount of oxygen incorporated into the aluminum, and clearly indicates that the oxygen 
penetrates several monolayers into the aluminum bulk. 

The solid black squares in Figure 12 represent a second run made three days after 
the first run saturated the surface with oxygen. The sample was not cleaned with an Ar + 
beam before the second run. This second run is characterized by a mass gain per incident 
ion, with an initial value that is only about three-fifths the initial value of the gain 
measured during the first run. We conclude that much of the oxygen incorporated into 
the aluminum during the first run did not remain near the surface. The fact that the 
second run was able to incorporate 8.9 nanograms of oxygen into a sample which had 
three days before been saturated with 13.8 nanograms of oxygen reveals that -60% of the 
originally incorporated oxygen had diffused out of the near-surface volume. We are not 
able to distinguish between diffusion of the excess oxygen deeper into the aluminum bulk 
and diffusion out from the surface into the vacuum. 

Figure 13 shows the mass change caused by two different C> 2 + beams of the 
same energy but of different intensities. The beam with the smaller current permitted 
accurate measurements of mass change in the region of the curve corresponding to small 
cumulative ion doses; and the beam with the larger current allowed accurate 
measurements for large doses. Note that in the region of intermediate doses, the two 
data plots coincide. This agreement is particularly noteworthy in view of the fact that the 
two runs were conducted on different aluminum samples. The agreement indicates that 
the time scale of the oxygen diffusion is either much larger than or much smaller than the 
time required to acquire the data for a single run. However, if the time scale were much 
smaller, then the mass change per incident ion would remain constant until oxygen 
diffused throughout the entire aluminum sample, contrary to the data. Therefore it is 
probable that the time scale for diffusion is long compared to the run time. This 
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conclusion is supported by the 60% "refilling" of the oxygenated sample by the second 
200-eV 0 + run three days after the first run. 

Figures 12 and 13 reveal another interesting aspect of the oxygenation process. 
When the aluminum contains very little oxygen (small dose), nearly every 200-eV 0+ 
ion that strikes the surface is embedded in the aluminum. In contrast, approximately 
sixty percent of the initial 200-eV Oj + ions is embedded into the clean aluminum. 

It is necessary for us to consider the effects of the ambient gases of the vacuum 
chamber. Since the pressure in the chamber increases by an order of magnitude when a 
beam is turned on, most of the ambient gas is from the ion source, primarily O 2 during 
the 0 + and C> 2 + bombardments. The influence of the ambient gas on the mass-loss 
measurement depends on the pressure, the sticking coefficient of the gas, and the current 
density of the ion beam (since this dictates the removal rate of adsorbed atoms and 
molecules). Below we develop a set of rate equations that describe the different 
processes which can either increase or decrease the mass of the sample. Note that 
molecular oxygen dissociates when the molecule adsorbs to the metal surface. 

In order to quantitatively estimate the effects of the ambient gas, we need to know 
the flux and the sticking coefficient as a function of coverage. For an aluminum sample 
exposed to O 2 gas at room temperature, data describing the coverage, 0, as a function of 
L (L = 10‘ 6 Torr sec) can be extracted from published experiments (Hayden et al, 1981 
and Yu et al., 1980 for example). It should be noted that, at this point, we are only 
discussing exposure of the aluminum surface to oxygen gas. Other than the normal 
thermal motion of the gas molecules, the gas is not accelerated toward the surface. 

This coverage of the sample by the adsorption of gases can be described by 
d0 1 

— = (1 - 6)— with initial condition 0(t ) I _ n = 0 (15) 

dt r ,-u 

when starting with a clean surface. For equation (15) and the following discussion, 1/x is 
the initial rate of adsorption for a clean surface. Stated differently, x is the incident flux 
multiplied by the sticking coefficient which obviously makes x dependent on the pressure 
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of the ambient gas. Solving equation (15), we obtain the coverage from the adsorbed gas 
as a function of time. 

-r 

6(t)=l-e~ 7 (16) 

-t 

where 0 = 1 corresponds to the saturated value. In general, 6(t) = 0 O (1 - e T ) where 0o is 
the portion of a monolayers coverage. A value for t can be obtained experimentally or 
from the literature. The information regarding the sticking coefficient as a function of 
coverage is included in equation (15). 

If in addition to exposure to the ambient gas the sample is bombarded with an ion 
beam, equation (15) should contain a term that characterizes the possible removal of the 
adsorbed gases by the ion beam 


dd_ 

dt 


= ( 1 - 0 )-- y 1 3 


T 


Nq' q 


e 


(17) 


where Yj is the sputtering yield for the adsorbed atoms when 0 = 1, J is the current 
density, and No' is the surface density of adsorbed atoms when 0=1. If we define No as 
the number of adsorbed atoms per unit area when 0=1, and for brevity define 


y ; 


N 0 q 


OB) 


Equation (17) can then be expressed as 

df) 1 

-f = d -e)--Kje. (19) 

dt T 

where I is the current that strikes a unit area of the sample. The general solution to 
equation (19) is 

;/ V 

6(t) = C-e ,i/ ' +Kan +- r /± — (20) 

1 A +K o 1 


where C can be determined from the initial conditions: 1) starting with a clean surface, 
6(0) = 0, and 


( 21 ) 
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or, 2) starting with a covered surface, 6(0) = 1, and 

6(t) = 77 7^ 77 k/- e'° A+Ka n +!/ T \. (22) 

(y T +K 0 I) L /T i 

It is clear from equations (21) and (22) that higher beam currents will decrease the 
time until an equilibrium exists between the adsorption of more gas and the removal of 
adsorbed atoms. The coverage when this occurs will be referred to as the equilibrium 
coverage, 6 eq , and is independent of the initial conditions. 




" (%+*<,!) 


Furthermore, if ff 0 I » A , 6gq is negligible. 

Of course we are primarily interested in the coverage that arises from the ambient 
gas because it changes the mass of the sample. The mass change of the sample as a 
function of current when the sample is being bombarded by an ion beam is described by 


dm dN 
= m n 


m AI • y 0 -(1 — 0 ) — 

4 


where dm/dt is the mass change of the sample occurring over the beam spot, mo is the 
mass of an adsorbed oxygen atom, m^i is the mass of an aluminum atom, and Yo is the 
measured sputtering yield when 6=0 (i.e. no coverage). Formally one should include 


the sputtering of oxygen atoms that are embedded in the volume by the primary beam, 
but one can never avoid the process of filling the volume and sputtering the oxygen. 
Therefore we note that our measured sputtering yield values are the effective sputtering 
yields of aluminum when the process of sputtering embedded oxygen is included in the 
measurement. In equation (24), N is the number of oxygen atoms adsorbed on the surface 
within the beam spot area, and can be expressed as a function of the coverage 

N(t) = N 0 6(t). (25) 


The first term in equation (24) describes the mass change associated with the 
adsorption of oxygen and the subsequent sputtering of these adsorbed atoms. The first 
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term can give a net positive contribution depending on the initial conditions, but is always 
zero when the equilibrium condition is reached. The second term describes the mass loss 
associated with the removal of aluminum, and is therefore always negative (because of 
the minus sign) for incident beam energies greater than the sputtering threshold. The 
second term is actually only an approximation since it suggests that the sputtering 
efficiency is proportional to the uncovered surface. Even for a covered surface, some 
sputtering of the substrate will occur, but if our experimental conditions lead to a small 
6 e q the approximation is quite good. 

We are now ready to substitute the coverage equation (either equation (21) or (22) 
depending on the initial conditions) into equation (24) to obtain: 


dm N n - + *o' 

— = — - — — • e L j — < 
dt T 


r 

1 

m A, • V 7 


r i ,i 



— +K 0 I 


... 

Lt J 

_ 




(26) 


for the case when the surface is initially clean, 6(0) = 0 


or 


dm 

hi 


= ~ m o N o 


K 0 I-e 



r 


> 

m Al ->W 



'i i 

’ 

q- 

~+K 0 I 



Lt J 

j 



(27) 


if the surface was initially covered, 6(0) = 1. 

To obtain the measured sputtering yield ( Y meas ), we only need to normalize dm/dt 
by the mass of an aluminum atom and by the number of incident ions. 


-dm 

Y (t) = -dt . 

± meas V 1 / J 

Q 

Both equation (26) and (27) reach the same equilibrium value 


(28) 


i.e. 


= V 0 - 0(OL J = VO- oj 


(29) 
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The two cases, however, approach the final value from opposite directions. Figure 14, a 
plot of equation (26), shows that the equilibrium value is approached from above when 
the surface is clean at the beginning of the run. 



Figure 14. Calculated mass change per incident ion as a function of bombardment time 
for different currents where the mass change is in units of the incident 
oxygen ion. 


Figure 15, a plot of equation (27), shows that the equilibrium sputtering value is 
approached from lower values when the surface is initially covered with adsorbed atoms. 
Irrespective of the initial conditions, the final measured sputtering yield at the 
equilibrium condition reaches the same value expressed in equation (29). We have an 
exponential behavior with characteristic time, teff, 


1 f l 

— = I —+K n I 


(30) 


If we substitute the value for 6 eq from equation (23) into equation (29), we obtain 

f 1 \ 


C„ = n {i-e H )=r, 


1 - 


1 


~+K 0 I 

V T ) 


= Y, 


K 0 I 


0 1 


(31) 


+ K 0 I 
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Figure 15. Calculated mass change per incident ion as function of bombardment time 
for different currents where the mass change is in units of the incident 
oxygen ion. 

which reveals that when K 0 I » — the final measured equilibrium sputtering yield is 

T 

approximately equal to the sputtering yield of the aluminum without the effects of the 
ambient gas. 

We are now in a position to estimate the magnitude of the effects of the ambient 
gas on our measurements using worst-case parameters from our experiment. 

f 2 particles Y 1Q -9 C \ 

Kq I = —1— = 1 iSS-A. secj — 1-sec-' (32) 

N « q (4.7 1 <)" particle. f 1 . 6 1 (T 19 — j 3800 


if the current were 1 nA. This means that a current of 1 nA removes 1/3800 of a 


monolayer per second. From the work of Hayden et al. (1981), we can ascertain that 


T 


2000 


-sec 


-i 


(33) 


for oxygen on aluminum. This indicates that even if our currents were as low as 7 nA 
(the smallest current we used), then the condition K 0 I » ~ is met and Y e f ea , » Y 0 


We can also calculate the time needed to effectively reach this equilibrium condition. 



38 


3t =3 (-+*•„/] =s( + - 

ff It ) v 2000 sec 3800sec 


This time is much shorter than was experimentally observed before the measured 
mass change per incident ion reached a constant value which signaled that the 
equilibrium condition had been achieved. We have shown that the adsorbed gas atoms do 
not appreciably affect our mass loss measurements after equilibrium conditions are 
reached, but the above description does not explain the very long time (i.e. large dose) 
that was necessary before reaching the equilibrium condition. To explain this ,we need to 
include the process involving the dose dependent mass gain associated with the 
implantation of projectile ions into the aluminum. The aluminum sample obviously 
incorporates oxygen into the surface and near-surface bulk since we observe a fairly large 
mass increase of the sample before we can measure any net erosion. The aluminum 
sample can incorporate oxygen most efficiently when it contains the least amount of 
oxygen. We will now show that the implantation of projectile atoms into the bulk is 
responsible for the long-term effects that we measure. 


To correct our rate equation for this process, we need to add the term 


m °' v ° li 


to the right side of equation (24) where Vo is the number of oxygen atoms implanted into 
the volume when the equilibrium condition is met (i.e. when volume is saturated), and V 
is the degree to which the volume is filled. The form of the equation that describes the 
rate that V changes with time is similar to the form of the coverage equation 


^ = d ~V)~ (36) 

dt V 0 q 

where y is a coefficient of entrapment that represents the probability that the oxygen ions 
will remain trapped in the near-surface bulk when V = 0. 

If the surface is initially clean, then V(r)| i=0 = 0 and the solution to equation (36) 


has the form 
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We then have 


V(t) 


dV y 1 

m o V o~r = m a V o ~ — e 
dt V 0 q 
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-yt i 

_ hill 

_V 0 q 
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(37) 


(38) 


for the volume term which is always greater than or equal to zero. Since this term 
reaches zero after the equilibrium condition is achieved, it does not affect our measured 
equilibrium sputtering yield. 

The mass change of the sample as a function of time can now be expressed if 
equation (38) is substituted into equation (26) if 6(0) = 0 or equation (27) if 6(0) = 1 i.e. 


dm m 0 N 0 
dt T 
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for the case when the surface is initially clean, 6(0) = 0 or 
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(40) 


for the case when the surface is initially covered, 6(0) = 1. A plot of equation (39) is 
shown in Figure 16 which agrees closely with the behavior seen in Figure 13 where a 
200-eV 02 + beam was incident on an aluminum sample that had been cleaned with 
argon. The bombardment time observed in our experiment before the mass change per 
incident ion decreased to zero is remarkably close to the value shown by the calculated 
curve. 

Our predictions are also quite close to the experimental observations for a surface 
that is initially coated before bombardment. A plot of equation (40) is shown in Figure 
17 which agrees closely with the behavior seen in Figure 12 where a 10 nA, 200-eV 0 + 
beam was incident on an aluminum sample which had been exposed to the ambient gas 
for a long enough period to at least partially cover the surface. Note that the rise that 
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occurs in the beginning of several of the oxygenation runs is also present in the calculated 
plot for the surfaces that start with at least a partial coverage. 



Figure 16. Calculated mass change per incident ion for the 75 nA, 200 eV 02 + run as 
a function of bombardment time if the volume effect is included. 


We can determine a value for Vo and y from our data, and then use a procedure 
similar to the one used for equation (24) to determine the time necessary to reach the 
equilibrium condition for the volume process, t Veg . We will show the calculation for the 
high current 200-eV C> 2 + run since we had enough current for that run to easily determine 
the time needed to reach the equilibrium condition. 


/ = 3t =3 

l Veq Vejf J 


u 

v 0 f 


v 


= 3 


(0.6) (75 - 10 ) 


(5.2- 10’ 4 )( 1.610 ,v ) 


= 1.5 hr. 


(41) 


This is nearly equal to the elapsed time observed during the experiment before the 


equilibrium condition was reached. A similar calculation for the low current (10 nA) 
200-eV 0+ run predicts that the equilibrium time should be approximately 7 hours. This 
was also very close to the time observed in that experiment. 
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Figure 17. Calculated mass change per incident ion for the 10 nA, 200 eV 0 + run as a 
function of bombardment time if the volume effect is included. 

Sputtering Yield 

For the aluminum sample, the Ar + sputtering yields are consistently larger than 
those for oxygen or nitrogen ions for a given energy. For purely momentum-transfer 
sputtering, one would expect N 2 + to give the highest sputtering yield because of the 
target-projectile mass ratio. However, other factors such as electronic interactions and 
chemical effects can influence the total sputtering yield. When the projectile species is 
reactive with the sample material, the chemical effects can either increase or decrease the 
effects of pure momentum-transfer sputtering. In chemical sputtering, energy of the 
projectile is first transferred to an electronic excitation of the substrate which can then 
impart energy to substrate particles causing their ejection from the surface. When 
chemical sputtering becomes an important contribution, the sputtering yield will be 
strongly dependent on the temperature of the sample. An increased sputtering yield is 
normally observed only at low energies where the effects of chemical sputtering can 
become more important. However if the reactivity of the bombarding species drastically 
increases the amount of the incident beam that remains embedded in the surface of the 
sample, the net sputtering yield normally decreases. For low energies, the decrease in 
sputtering yield can be attributed to two effects: 1) the altered chemical composition of 
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the surface layers often leading to larger effective surface binding energies, and 2) the 
fact that the incident particles must sputter these incorporated projectile atoms in addition 
to sputtering atoms from the sample material (Roth, 1983). 

A comparison of the measured sputtering yield of 2-keV Ar + on a clean aluminum 
surface (5.4xl0 -23 grams per incident ion) versus the sputtering yield of 2-keV Ar + on an 
oxygenated aluminum surface (2.2x1 O’ 23 grams per incident ion) reveals that the 
chemical composition of the surface layers is indeed an important factor. The oxygen 
that is mixed with the aluminum in the surface layers likely increases the effective 
binding energy of the surface which leads to a smaller sputtering yield. For the 2-keV 
case, the altered surface layers may also affect the collisional cascade which could result 
in a decrease in the sputtering yield. 

A decrease in the effective binding energy of the sample surface caused by the 
altered chemical composition of the surface layers may explain why the sputtering yields 
of 300- and 400-eV N2 + on the aluminum oxide sample are larger than the corresponding 
sputtering yields of C> 2 + . Momentum-transfer sputtering for oxygen is favored over that 
for nitrogen because of the match between the projectile and target masses. Despite these 
considerations, the sputtering yields for nitrogen were larger than those for oxygen, 
indicating that the altered chemical composition of the surface layers is important in this 
instance. 

The sputtering yield results of 2-keV Ar* incident on a clean aluminum surface, 
an oxygenated aluminum surface, and on an aluminum oxide surface are summarized in 
Table 6. The measured sputtering yield is largest for the pure aluminum case. It is of 
some surprise that the sputtering yield of argon on aluminum is approximately 40% 
smaller than it is for argon on aluminum oxide. A similar trend is seen from a 
comparison of the sputtering yields for 200-eV C> 2 + ions on an oxygenated aluminum 
surface to the yields for 200-eV C> 2 + ions on an aluminum oxide. The measured 
"equilibrium" sputtering yield of 200-eV C> 2 + on the oxygenated aluminum sample is 
again roughly 40% lower than the measured sputtering yield of 200-eV C> 2 + on the 



43 


aluminum oxide sample. This seems to suggest that the excess oxygen present in the 
oxygenated aluminum surface serves either to increase the effective binding energy on 
the surface or reduces the collisional cascade thus decreasing the sputtering yield. 


Table 6 . Comparison of sputtering yields from 2 keV Ar + on similar surfaces. 


Projectile Energy 

Projectile 

Sample 

Y 

(grams lost per incident ion) 

2 keV 

Ar+ 

A1 

5.4x1 0 - 23 

2 keV 

Ar+ 

oxygenated A1 

2 . 2 xl 0- 23 

2 keV 

Ar+ 

AI 2 O 3 

3.7x10-23 




SUMMARY AND CONCLUSIONS 


This report describes experiments designed to investigate the interactions of low- 
energy reactive ions with surfaces and measure the sputtering yields from low-energy ion 
beams (oxygen and nitrogen). We constructed an experimental apparatus to produce a 
low-energy ionic oxygen or nitrogen beam which we used to bombard aluminum and 
aluminum oxide samples, measuring the mass change of the sample per incident oxygen 
ion. A novel technique for increasing the sensitivity of a quartz-crystal microbalance by 
a factor of six allowed us to measure the very small mass changes induced when low- 
current, low-energy ion beams were incident on the thin film samples 

Certain aspects of the experiment were of particular interest. For instance, when 
200-eV C>2 + or 0 + ions are incident on a clean aluminum surface, the mass of the sample 
increases until the surface has been exposed to a large dose of oxygen. For small 0 + dose 
(when the surface and near-surface bulk contain little oxygen), nearly every 0 + ion 
becomes embedded in the aluminum. For small C>2 + dose, only 6 out of every 10 C>2 + 
ions that strike the surface are embedded. This oxygenation process alters the chemical 
composition of the surface and the near-surface bulk. 

One of the most striking aspects of the oxygenation process concerns the total 
amount of oxygen that the clean aluminum surface and near-surface bulk can hold. Our 
experiments revealed that the total mass increase of the aluminum is the same for 
bombardment by 200-eV 02 + or 0 + ions. We calculate that the mass increase per square 
centimeter of surface area is equivalent to a minimum of 1 1 oxygen atoms for every 
aluminum atom on the surface. A 200-eV 0 + run on an aluminum sample that had been 
saturated with oxygen three days earlier revealed that approximately 60% of the 
embedded oxygen diffuses out of the saturated volume after several days. Most of the 
remaining oxygen is probably in the form of AI2O3. 
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Rate equations were derived that describe the contributions of the five competing 
processes that change the mass of the aluminum sample during these experiments: 1) 
sputtering of the aluminum, 2) incorporation of oxygen from the ion beam into the 
aluminum, 3) sputtering of incorporated oxygen that came from the beam, 4) 
incorporation of oxygen from the background gas in the sample chamber, and 5) 
sputtering of the adsorbed background gas atoms. With these rate equations, we are able 
to calculate the maximum amount of incorporated oxygen and the time (for a given beam 
current) before the mass change per incident ion becomes a constant negative value 
(signifying oxygen saturation). These calculated values agree remarkably well with the 
experimental results. 

Sputtering yield measurements were made for low-energy 0 + , C> 2 + , N + , and N2 + 
on copper, aluminum, and aluminum oxide samples, and 0 + and C> 2 + on silicon oxide 
samples. The sputtering yields of copper under oxygen bombardment exhibited time- 
dependent behavior similar to that for aluminum, which is consistent with the fact that 
both form stable oxides. The sputtering yields for oxygen bombardment of aluminum 
decrease as the oxide layers form, giving the aluminum a protective coating that is more 
resistant to erosion. The sputtering yields measured for silicon oxide were very small, 
even for 2-keV oxygen beams, indicating that silicon oxide is very stable in a reactive 
oxygen environment. 

These results will help quantify the rate that orbiting vehicles will erode when 
exposed to the ambient space environment. The results also suggest that the altered 
chemical composition of the surface that occurs plays an important role in the sputtering 
process, particularly for oxygen bombardment of aluminum and nitrogen bombardment 

of aluminum oxide. 
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